Breathing mode of rapidly rotating Bose-Einstein condensates 
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We show that the breathing mode of a rapidly-rotating, harmonically-trapped Bose-Einstein con- 
densate may be described by a generalized lowest Landau level (LLL) wave function, in which the 
oscillator length is treated as a variable. Using this wave function in a variational Lagrangian for- 
malism, we show that the frequency of the breathing mode for a two-dimensional cloud is 2uj±, 
where ui±_ is the trap frequency. We also study large-amplitude oscillations and confirm that the 
above result is not limited to linear oscillations. The resulting mode frequency can be understood 
in terms of orbits of a single particle in a harmonic trap. The mode frequency is also calculated 
for a cloud in three dimensions and the result for the axial breathing mode frequency agrees with 
recent experimental data in the rapid rotation regime. 

PACS numbers: 03.75.Kk, 05.30.Jp, 67.40.Vs, 47.37.-|-q 



I. INTRODUCTION 

The creation of vortices and vortex lattices in Bose- 
Einstein condensates of cold atomic gases Q, H, |^ has 
opened up the study of vortex lattices in a regime which 
cannot be studied in liquid helium 4. Due to the dilutc- 
ness of cold atomic gases, it is possible to realize situa- 
tions in which the vortex core size is comparable to the 
vortex spacing 4] . 

In a seminal work, Ho pointed out that the Hamilto- 
nian for a rotating cloud in a harmonic trap has the same 
form as that for charged particles in a uniform magnetic 
field. Based on the similarity between these two systems, 
he argued that when the rotation frequency Q. is close to 
the transverse frequency u>^_ of the trap, almost all parti- 
cles would condense into the lowest Landau level (LLL) 
of the Coriolis force. Stimulated by this insight, exten- 
sive experimental studies have been performed by the 
JILA group (see, e.g., Coddington et al. j^), who have 
achieved angular velocities in excess of 0.99wj_ at which 
the cloud contains a vortex lattice with several hundred 
vortices 0. 

The current frontier of the experiments is entering the 
mean-field LLL regime in which fiVt is comparable to or 
larger than the interaction energy, gn, leading to a system 
whose wave function is dominated by the LLL compo- 
nent. Here n is the particle density and g = ATrh^as/m is 
the two-body interaction strength, m is the particle mass, 
and flg is the s-wave scattering length. In this regime, the 
Gross-Pitaevskii equation is still applicable because the 
number of particles N is much larger than that of vor- 
tices A'v. As shown by ourselves Ig and confirmed by 
Cooper et al. ;9j (see also Ref. T(T|), a single-particle 
wave function of the LLL form is a good approximation 
to the ground state in the rotating frame, i.e., the yrast 
state, of rapidly rotating condensates. Even the subtle 
connection between the vortex lattice distortion and the 
density profile can be described by the simple LLL ap- 
proximation: a small distortion of the vortex results in 



a drastic change of the density profile from Gaussian to 
the Thomas-Fermi parabola. 

However, limitations of the LLL wave function become 
apparent when we consider excited states. The oscilla- 
tions of the cloud radius (a so-called breathing mode) in 
the direction of the rotation axis have been measured in 
the mean-field LLL regime 7j, and those in the trans- 
verse direction are expected to be measured in the fu- 
ture experiments. However, the LLL wave function does 
not have the fiexibility to describe these simple collec- 
tive modes. In the present work, we consider which de- 
grees of freedom in the LLL wave function are relevant 
to describe the breathing mode and provide a simple im- 
proved treatment, which remedies the limitations of the 
LLL wave function. 

In the beginning of the next section, we give the basic 
formalism of the present analysis: the variational La- 
grangian approach Til. Using this framework, we then 
discuss the breathing mode of the rapidly rotating two 
dimensional Bose-Einstein condensates. In Section IIIII 
large-amplitude oscillations are discussed. In Section llVl 
we provide a simple physical understanding of the result. 
Summary and conclusion are given in Section The 
mode frequency in three dimensions is calculated in Ap- 
pendix^ 



II. GENERAL FORMALISM 

Let us first see why the standard LLL wave function 
cannot describe the breathing mode. For definiteness, we 
consider purely two-dimensional motion in the plane per- 
pendicular to the rotation axis. We write the condensate 
wave function ij^ in the form ^ — \/N(f), where (p is nor- 
malized such that integrated over all space is unity. 
The LLL wave function can be written as 
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where the rotation axis is along the 2:-axis, C = x -\r iy, 
Q are the vortex positions measured in the laboratory 
frame, A'v is the number of vortices, r = (x, y), and d±^ = 
^Jhjmtjj\_ is the transverse oscillator length. 

Employing Eq. we find the angular momentum per 
particle to be given by 



rf2 



(2) 

where (|(/>lllP) is the density profile smoothed over an 
area of linear size large compared to the vortex separa- 
tion but smaller than the radial extent of the cloud. For 
a uniform density of vortices, the density profile has a 
Gaussian form. 
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where a is the width parameter given by the vortex den- 
sity riv as riv = l/(7rd^) — l/(7rcr^) In this case one 
finds 



(4) 



A. Extended Lowest Landau Level Wave Function 

To describe the breathing mode in the rapidly rotat- 
ing limit, we introduce an extended lowest Landau level 
wave function, which keeps the LLL form but allows the 
oscillator length of the trap (i.e., the trap frequency) to 
be adjustable: 



0cx(r) = Aox ]^(C - 0) exp 



1 \ 



(7) 



The dynamical variable A describes the variation of the 
particle density. The wave function (Q, like the original 
LLL wave function has a Gaussian smoothed density 
profile if the lattice is uniform, but with a modified width 
parameter given by 



(8) 
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The other dynamical variable /3 generates a radial ve- 
locity field Vr which causes a homologous change of the 
density profile |l2l | and the corresponding velocity is 



In the real situation, interactions distort the lattice and 
the density profile becomes a Thomas-Fermi parabola |^ 



Vr 



h d 
m dr 
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m dp. 



(9) 



{\<PiLV)^^{r) = m{i 
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with 1^(0) — 2/tiR\, where R±^ is the cloud radius. In 
this case the angular momentum is 



Ml 



(6) 



Here the superscripts "G" and "TF" denote that the den- 
sity profile is Gaussian and Thomas-Fermi, respectively. 

Equations l@J and © show that, for the LLL wave 
function, the cloud width parameter a or the cloud radius 
is fixed if the angular momentum is constant, i.e., the 
breathing mode cannot be described by the standard LLL 
wave function. 



In the present analysis, we use the variational La- 
grangian formalism 11]. The Lagrangian per particle 
consists of the time-dependent part T and the energy 
functional E: 



with 



Eint [00 
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where V, and Ei^t are the kinetic, potential, and inter- 
action energies. Each term in the Lagrangian functional 
for the extended LLL wave function can be written as 



dt 



dt 
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and 
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where (72D is the effective interaction parameter in two 
dimensions. If the system is uniform in the z-direction, 
32D = Ng/Z, where Z is the axial extent of the cloud, 
while if the system in the z-direction is in the ground 
state of a particle in a harmonic potential of frequency 
LUz, then (72D = Ng/ {dzV^Tr), where d^ = y/hlmiOz is 
the oscillator length in the z-direction [l^- The angular 
momentum per particle for (/)ex is 



• X (-iW^)]^ 
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Let us now consider the first term of Eq. H12I) , which 
we denote by Ti . The velocity of the vortex in the lab- 
oratory frame has the azimuthal component v^^i ~ r^rj 
due to the bulk rotation and the radial component Vr^i — 
6R ujTi /R, where R is the scale of the radial extent of the 
cloud, SR and to are the amplitude and the frequency of 
the breathing oscillation. Employing the averaged vortex 
approximation as we will do in the later discussion, only 
the azimuthal velocity gives a nonzero contribution be- 
cause the coarse-grained density profile is axisymmetric: 



2 (r - r^) X 4,,v^,, 
d r ■ --—^v{r) , 



(17) 



where (f)^ is the unit vector of the azimuthal direction 
measured from vortex i. We can readily see that the in- 
tegral of Ti vanishes when the density is uniform. For 
the inhomogeneous cloud, the dominant contribution of 
Ti comes from vortices at ^ _R; vortices at <C i? 
give only small contribution because the cancellation of 
the cross product between different r is efficient [l^ . 
The integral of Ti for each vortex at ^ _R yields 
~ / d^r Q,v{r) — Q. because |r — r^j ~ i? and v^_i ^ RQ,. 
The number of vortices in the surface region is propor- 
tional to n^lR ^ Nvi/R, where tt£'^ gives the area per 
vortex. Thus, for a cloud with a large vortex lattice, Ti 
scales as 



Ti cx nnN^ 



R 



(18) 



which becomes negligible compared to the second term 
of Eq. (Hai T2 ~ hiuN^ if i? > ^, i.e., N^^ > 1 [we will 
see later that uj ~ r2(~ (^±)]- In the later discussion, we 
assume that the vortex lattice is large enough that the 
effect of Ti may be neglected. 



B. Gaussian Profile 

Let us first consider a cloud with a uniform vortex 
lattice. Using the averaged vortex approximation, we 
replace |0cxP in Eqs. (|12|l - IjlSI) by the smoothed density 
profile of a cloud ||2J) with a uniform vortex lattice. We 
then obtain 



(19) 

-f3')^, (20) 
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where b = (|?!'ox|'*)/(|0exP)^ is the Abrikosov parameter, 
which is comparable to unity. 
Similarly, Eq. l(TB|l reads 
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Due to the angular momentum conservation, the dynam- 
ical variables a and A are not independent by Eq. H23|l . 
Writing the constant value of Lz as Iz, \ is expressed as 



1 

A2 



dl 
^2 
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(24) 



Using Eq. (|24|l to eliminate A from the Lagrangian, we 
obtain 
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Thus the Euler-Lagrange equations for the dynamical 
variables [3 and a 



dt\ dp J dp 



dt \ da J da 
lead to the following equations, 
a = w^/3(T 

and 



= , 
= , 



+ r + 1 
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respectively; where k = mbg2D / i'^Trh^) is the dimension- 
less interaction strength. Combining the above two equa- 
tions, we obtain an equation of motion of the cloud width 
parameter a: 



(30) 



According to Eq. 1)30(1 . the cloud width parameter ctq in 
the equilibrium state is given as 



diif 



(31) 



We note that this expression coincides with the width 
parameter of the minimum energy state obtained by a 
variational calculation in the rotating frame. Now we 
suppose a oscillates around ctq ^-s 



a{t) fTo + Sa{t) 



(32) 



and its deviation Sa{t) is small. Linearizing Eq. H30|) with 
respect to Sa yields 



Sa{t) + {2uj^_)'^Sa{t) = 



(33) 



Thus the breathing mode frequency to in the LLL regime 
for a Gaussian profile is 



ui = 2ijj\ 



This result agrees with what one obtains from calcula- 
tions based on hydrodynamics if one takes the limit of a 
trapping potential that is independent of z .15> .16.J . In 
the case of Ref. |0, the polytropic index must be put 
equal to the value two appropriate for a dilute Bose gas. 



It is notable that the interaction energy scales in the 
same way as the rotational kinetic energy (the P term) 
in the Lagrangian 1(25(1 and in the equation of motion 
Thus the interaction energy enters only in the 
combination P + k, so it does not affect the frequency 
of the breathing oscillation This is a remarkable 

feature of the two-dimensional system ^3 ■ 

In the case of a three-dimensional rotating cloud 
trapped in a potential V — mLo±{x^ -I- y^)/2 -I- mijJzZ^ jl^ 
the interaction energy is given as E\^x_ ^ l/{a'^Rz), 
where is the radius (width parameter) of the cloud 
in the z-direction when the density profile of this direc- 
tion is the Thomas- Fermi parabola (Gaussian). Thus, 
for a non-zero interaction strength, the two-dimensional 
monopole oscillation in the xy-plane couples with the 
one-dimensional monopole oscillation along the z-axis 
whose frequency is v^Wz (in the case where the zero- 
point energy of the z-direction is negligible compared to 
the interaction energy) and the resulting mode frequency 
is modified from 2uo± (see Appendix 



C. Thomas- Fermi Profile 

As in the previous section, we adopt the averaged vor- 
tex approximation. For the inverted parabolic density 
profile of Eq. © for a distorted vortex lattice, Eqs. ifT^ 
- ((T3|) can be written as 
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cx J 
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huj± R\ 
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(37) 
(38) 



Using the angular momentum conservation: 

1 I , (39) 



(34) A can be expressed as 



1-3^ ('^ + l\='i^l 



(40) 



Thus the Lagrangian can be written with only i?^ and 
/?: 



5 



h f3Rl 
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The Euler-Lagrange equation for (3 leads to 
and that for i? 1 to 



Combining Eqs. (|42|l and 143() . we obtain 
1 

+ 



- 9 P 

d1 



Landau levels, i.e., gn ^ ?ia;^ or Ngb/Z ^ R'^/d\, 
where R denotes the cloud width or radius. To keep 

(42) the two-dimensional character of the system, we require 
hujz S> mR^; otherwise, oscillations in the z-direction are 
excited. One should also note the criterion for the first 
term Ti in the time-dependent part of the Lagrangian 

(43) l(T^ to be neghgible compared to the second term T2. 
If the amplitude is large, the cloud can shrink signif- 
icantly as the angular velocity (and Ti) can become 
large. Feynman's relation gives the angular velocity of 
the cloud as ~ h/{m£'^) ~ fiN^ / {mR^)\ then we have 



K = . (44) Ti - nnN. 



The cloud radius i?^,o in the equilibrium state is thus 
given as 



(45) 



which coincides with the cloud radius obtained by mini- 
mizing the energy in the rotating frame. Again we notice 
that the interaction energy scales in the same way as 
the rotational kinetic energy and it does not affect the 
breathing mode frequency |l8j| (see the Appendix for the 
three-dimensional case). 

A linearized equation of motion for a small oscillation 
oiR^it) = R±^o + SR±{t) is 



SR±it) + i2u±f SR±{t) = 



(46) 



According to Eq. H46() , the breathing mode frequency u in 
the LLL regime for a cloud with a Thomas-Fermi profile 
is also 



± 



(47) 



which coincides with results in the slow rotation regime 
[TEL fl^ . It is also noted that this result is in accord with 
measurements in the slowly rotating regime. 



III. NON-LINEAR OSCILLATIONS 

In the present section, we consider large-amplitude os- 
cillations. However, we should mention that the sev- 
eral assumptions in the preceding discussion can break 
down in the non-linear regime even if the mean-field 
Gross-Pitaevskii theory is still a good approximation, 
and thus the validity of the result in the present sec- 
tion is limited. First of all, for the extended LLL wave 
function to be valid, the density should be always low 
enough to ensure that the interaction energy is much 
smaller than the energy gap between the LLL and higher 



fi^ Ny^"^ / {mR^). Thus the criterion for 

1/4, 



Ti < r2 - huj^N^ leads to i? > K' d^. 

The equations of motion H30|l and (|44() for both the 
Gaussian and Thomas-Fermi profiles can be written as 
the following general form: 



'dp 



d 
dX 



2 



2A2 



0, 



(48) 



where X denotes a or Xq corresponds to ctq or 

R±,o given by Eqs. (|31|l and H45|) . respectively. Here 
we measure the length and time in units of d± and 
(^J^- Equation H48(l has the same form as the Newton's 
equation of motion for a particle moving in a potential 
V{X) EE X^/2 + A4/2A2 (see Fig. P). This potential 
shows the restoring force caused by the first term when 
the cloud expands and the strong centrifugal repulsion 
by the second term when the cloud contracts. 

Suppose that, at i — 0, the cloud has its equilib- 
rium size X = Xq but with velocity ^(0) — Xq to 
excite a breathing oscillation. Thus the initial kinetic 
and potential energies are Xq/2 and V{Xo) = Xq, re- 
spectively. The oscillation enters the non-linear regime 
when Xq/2 > Xq because the potential is no longer ap- 
proximated by the harmonic one. We thus introduce the 
"non-linearity parameter" x as 



X 



ii. 

2Xi ' 



(49) 



X ^ 1 in the linear regime and x ^ 1 in the non-linear 
regime. 

Multiplying Eq. by X and writing as XX = 

{l/2){dX'^/dt) and X{dVldX) = dV/dt, we can inte- 
grate this equation analytically and finally obtain 



dX 



^2{E-V{X)} 



= =F — arctan 



E-X^ 



X^2{E-V{X)} 



(50) 
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FIG. 1: Effective potential V{X) = + Xtj2X'^ of the 

motion of the radial extent of the cloud X. 



or 



tan [T2(i - to)] = 



E-X' 



X^2{E-V{X)} 



(51) 



Xq are the 



where to and E = X^/2 + V{Xo) ^ X^/2 
integration constants. We can see from Eq. (|51(l that the 
period r of the oscillation is r = tt (in units of ujj^) and 
thus the mode frequency is 2t:/t — 2t: / (ttuj^^) = 2uj±_ 
even in the non-linear regime. 

In Fig. [3 we show the time evolution of the width 
parameter a and the radius R± of the cloud for some ini- 
tial conditions. In this plot, we set I = 100 and k = 100, 
which are values appropriate for recent experiments 0. 
The radial extent of the equilibrium state corresponding 
to these values of I and k is ctq — 10.02di^ for the Gaus- 
sian profile and R±fl — 17.36d± for the Thomas-Fermi 
one. Figure El^a) shows linear oscillations with x < 1 
and the amplitudes of X are almost symmetric above 
and below Xq in the both cases. Figure |2Ib) shows a 
case in the non-linear regime with x = 0(10). We ob- 
serve that the oscillations are asymmetric and this fea- 
ture is more prominent in the Gaussian case than in the 
Thomas- Fermi one. In Fig. |2Ic), we show strongly non- 
linear oscillations with x > 100 as a demonstration (how- 
ever, Ti is no longer negligible in this case). We can see 
that the mode frequency uj = 2lu± even in the strongly 
non-linear regime. 
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FIG. 2; (Color online) Time evolution of the width param- 
eter G and the radius R± of the cloud from the linear (a) to 
the strongly non-linear (c) oscillations. 



for the two dimensional rotating cloud in the simplest 
way by focusing on the orbit of a single particle. Sup- 
pose a particle of mass m moves in a circular orbit of 
radius r. The frequency of oscillations when the cir- 
cular particle orbit is perturbed can be interpreted as 
the breathing mode frequency of the cloud. Now we 
assume that each particle moves in accord with the 
cloud motion, i.e., a = r/X is constant, where X de- 
notes a or i?^ in the case of the Gaussian profile or 
the Thomas-Fermi one, respectively. In terms of the 
variable a, the density of the cloud can be written as 
n(r) = n(0) exp (— r^/cr^) = TVa^ exp (— a^)/(7rr^) for the 
Gaussian profile and n{r) = n{Q){l-r'^ / R\) = 2Na'^{l- 
a?)/{'Kr'^) for the Thomas-Fermi one. The interaction 
energy of a single particle is given by g2T>n(r) = 
where 7 = (7V(72Da^/7'') exp (— a^) in the Gaussian case 
and 7 = (2iV(72Da^/7i')(l — a^) in the Thomas-Fermi case. 
The energy of the particle is 



IV. PHYSICAL INTERPRETATION OF = 2uji_ 

In the preceding discussion, we have seen that, in the 
rapidly rotating regime, the frequency of the breath- 
ing mode of a two dimensional cloud is w = 2u;j_ as 
in the slow rotation regime [TH ITgI l20l | (and in 
the non-rotating case l21|). Furthermore, this result is 
not limited to small amplitude oscillations. We can un- 
derstand the above robustness of the result u — 2uj\ 



E : 
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2mr^ 



1_ 



(52) 



where v is the azimuthal velocity and Lz — rmv is the 
angular momentum of the particle, which is a constant of 
motion. The radius rp of the unperturbed orbit is given 

by dE/dr\ = -Ll/{mrl) + muj\rQ - 2'y/r^ = : 



{Ll + 2m7)i/2 



mu!± 



(53) 
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FIG. 3: Schematic picture of the energy of a particle rotating 
in a harmonic trap of the frequency u)±_ . Oscillation frequency 
of the radius of the orbit is 2ijj±_ . 



have seen that the modulation of the oscillator length in 
the original LLL wave function is an essential degree of 
freedom to describe the breathing mode. We have also 
studied large-amplitude oscillations; we have observed 
that the amplitude becomes asymmetric due to the non- 
linear effect, but the frequency is still 2lo^. Finally, we 
have provided a simple and physical understanding of the 
result Lo = 2ujj_ and its robustness for two-dimensional 
rotating clouds. 
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We then give a small perturbation Sr of r and obtain 
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L'i 



2mrQ 



1_ 
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1/2 , 



Sr' 



(54) 



which is shown schematically in Fig. |21 The first term 
gives the energy of the unperturbed motion; the second 
term shows that the frequency of the radial oscillation 
equals to 2uj±. Note that the interaction energy, which 
has the same r^^ dependence as the rotational kinetic 
energy, just gives a correction to the equilibrium state 
and it does not affect the dynamics. Thus the breathing 
mode frequency directly reflects the general property of 
the rotating single particle motion in two dimensions. 



V. SUMMARY AND CONCLUSION 

In this work, we have investigated the breathing os- 
cillation of a rapidly rotating two-dimensional Bose- 
Einstein condensate. Using the variational extended LLL 
wave function, which incorporates the change of the vor- 
tex density, we have shown that the mode frequency is 
Lu = 2uj±, as in the slow rotation regime. It would be 
valuable to confirm this prediction in future experiments 
on rapidly rotating Bose-Einstein condensates. There we 



APPENDIX A: BREATHING MODE IN THREE 
DIMENSIONS 

We consider the breathing mode in the case of three 
dimensions and see how the mode frequency is modified 
from 2lli±. In the following analysis, we assume that the 
radial density profile is the Thomas-Fermi parabola and 
the axial one is Gaussian. Thus the extended LLL wave 
function in this case is 



,(r) = Ae,^(C-C^) exp 



2 1 



1 \ r 



X exp 
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,2 1 



17 



2d? 



(Al) 



where cr^ is the width parameter of the axial density pro- 
file and 7 is the dynamical variable describing the veloc- 
ity field, which causes the homologous change of the axial 
density profile. The coarse-grained density profile can be 
written as 



KO) 1 



(A2) 



with i^(0) = 27r"3/2/^^2^^^)^ rp^g terms in the La- 
grangian per particle for the above wave function may 
be calculated to be 
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N, 



^ 2(r- Ti) X Ti 
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with (73D = Ng. The angular momentum per particle is which is the same as Eq. ()39|l . Thus the Lagrangian 

functional can be written as 



V3A2d2 



- 1 



1 



(A7) 



3/2 ,,2 -.^^1 ncj^ f 1 . , 



X2 



+ 1)^ + ^ 77772 + (7^ + 



2Z2 



3 X^Z 



, (A8) 



r 



where X = R±/dj_, Z — cr^/d^, and 
mbgsu 1 
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is the dimensionless interaction strength. 

The Euler-Lagrange equations for (3 and 7 lead to 



X = u^fiX , 
Z = u^_^Z , 



(AlO) 
(All) 



respectively. Using these equations, the equations of mo- 
tion for X and Z yield 
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The values of X and Z in the equilibrium state, 
Xq and Zq, are given by the following coupled equa- 
tions X^ - 9P - SksbZq^ = and Z^ - Z^^ - 
{SiK3d/5){uj±/uJz)Xq'^ — 0, and linearized equations for 
small deviations 6X and 6Z from Xq and Zq are 



1 A 7 

— X^jl -f AX'^SX + 8k3d = , (A14) 



and 



Lz'^5Z + 3Z',SZ+^-^ 



I6K3D uji_ 5X 



LUz Xq 



= . (A15) 



Now writing SX = Axe^'^* and SZ = Aze^'^*, and using 
the condition that Eqs. ljA14|l and ljA15|l have a non- 
trivial solution, we finally obtain 
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We note that the 1/Zq terms have come from the zero- 
point energy in the z-direction. From a variational cal- 
culation in the rotating frame, we find 9/^ = Xpr^Q/oj^^, 
where Qq is the angular velocity of the cloud in the equi- 
librium state for a given angular momentum. If one ne- 
glects the zero-point energy in the z-direction (23|, the 
above expression reduces to 

(A17) 

which is exactly the same as Eq. (45) of Ref. (for 
a polytropic index equal to two) and Eq. (10) of Ref. 
[l6j derived within the hydrodynamic theory in the slow- 
rotation regime. This agreement justifies our results in 
the present paper obtained with the extended LLL wave 
function even for the slow rotation regime. In the breath- 
ing mode, only the coarse-grained density and the aver- 
aged vortex density are relevant degrees of freedom when 
the number of vortices is large; the local vortex structure, 
which cannot be described by the extended LLL in the 
slow rotation regime, is irrelevant. Unlike the ordinary 
LLL wave function, the extended LLL wave function can 



describe the averaged vortex density correctly also in the 
slow rotation regime due to the extra degree of freedom 
A, which is the reason of the above agreement. (In the 
slow rotation regime, Ti <^ T2.) 

In the rapid rotation limit, where the 2-dependence of 
the wave function corresponds to the ground state of a 
particle in a harmonic potential, Zq ~ 1, and Eq. (|A16p 
leads to 



(A18) 

In the limit of CIq lu_i_, the two frequencies become 
Lu = 2uj±_ and uj = 2u)^_. The former value corresponds 
to the transverse breathing mode and the latter to the 
axial one. Unlike those for the hydrodynamic models 
0,^1, our calculations are therefore able to explain the 
experimentally obtained change in the axial breathing 
mode frequency from the value lo = V^lJz given by Eq. 
(|A17p to 2u!z when the interaction energy per particle 
falls below huz (see Ref. [H for details). 
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